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ABSTRACT 


A  case  study  of  the  Andreas  et  al.  (1984)  data  on  atmospheric  boundary 
layer  modification  in  the  marginal  ice  zone  is  made.  Our  model  is  a  two- 
dimensional,  multi-level,  linear  model  with  turbulence,  lateral  and  vertical 
advection,  and  radiation.  Good  agreement  between  observed  and  modeled 
temperature  cross-sections  is  obtained.  In  contrast  to  the  hypothesis  of 
Andreas  et  al.,  we  find  the  air  flow  is  stable  to  secondary  circulations. 
Cloud  top  longwave  cooling,  not  an  a ir- to-surf ace  heat  flux,  dominates  the 
cooling  of  the  boundary  layer.  The  accumulation  with  fetch  over  the  ice  of 
changes  in  the  surface  wind  field  are  shown  to  have  a  large  effect  on 
estimates  of  the  surface  wind  stress.  We  speculate  that  the  Andreas  et  al. 
estimates  of  the  drag  coefficient  over  the  compact  sea  ice  are  too  high. 
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Figure 


I.  A  schematic  diagram  of  a  two-dimensional  marginal  ice  zone.  A 

model  of  the  atmospheric  boundary  layer  including  simplified 
treatments  of  low  level  stratus  clouds  and  sea  ice  thermo¬ 
dynamics  is  used  to  model  the  modification  of  the  prescribed 
inflow.  Modeled  processes  include  turbulence,  lateral  and 
vertical  advection  and  radiation.  The  bulk  exchange  coef¬ 
ficients  for  momentum,  heat  and  moisture  are  functions  of  the 
ice  concentration. 

2a.  A  comparison  of  observed  (from  Andreas  et  al.,  1984;  solid)  and 
modeled  (dashed)  ice  edge  wind  profiles.  The  modeled  profile  is 
the  prescribed  inflow  of  the  model  (adapted  from  Andreas  et  al., 
1984) . 

2b.  A  comparison  of  observed  (from  Andreas  et  al.,  1984;  solid)  and 
modeled  (dashed)  ice  edge  potential  temperature  soundings.  The 
modeled  sounding  is  the  prescribed  inflow  of  the  model  (adapted 
from  Andreas  et  al.,  1984). 

2c.  The  temperature,  T,  dew  point,  TD,  and  liquid  water  mixing 

ratio,  r^,  soundings  used  to  prescribe  the  model  inflow. 

3.  Potential  temperature  soundings  at  the  ice  edge  and  at  fetches 

of  40,  80,  120  and  150  km  for  the  rough  simulation. 

4.  A  height-fetch  cross-section  of  potential  temperature  (°C; 

above)  and  the  initial  fraction  of  open  water  (below)  for  the 
rough  simulation.  The  stippled  regions  in  the  upper  figure 
indicate  liquid  water  mixing  ratios  greater  than  0.01  g/kg. 

5.  The  heat  budget  for  the  rate  of  change  of  potential  temperature 
of  the  air  column  at  a  fetch  of  150  km  in  the  rough  simulation. 
The  corresponding  potential  temperature  sounding  is  shown  in 
Figure  3.  Processes  include  turbulence,  total  radiation 
(longwave  plus  shortwave)  and  vertical  and  lateral  advection. 

6.  The  25  m  drag  coefficients,  CD,  and  bulk  exchange  coefficients 
for  sensible  heat,  C^,  for  the  rough  (solid),  smooth  (dashed), 
and  intermediate  (dotted)  simulations.  In  the  intermediate 
Simula  tion  CH  “  CD* 

7.  The  25  m  surface  wind  speed,  v  ,  and  surface  wind  stress, t  ,  of 
the  rough  (solid),  smooth  (dashed),  and  intermediate  (dotted) 
Simula  tions . 

8.  Same  as  Figure  4  but  for  the  intermediate  surface  simulation. 
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Introduc  tion 

The  marginal  ice  zone  (MIZ)  is  a  boundary  zone  between  the  polar  and 
temperate  parts  of  the  climate  system  (Un ters teiner  et  al.,  1983;  Johannessen 
et  al.,  1983).  An  important  aspect  of  the  MIZ  of  the  Bering  and  Weddell  Seas 
is  the  intense  rate  of  destruction  of  sea  ice  that  had  formed  in  the  coastal 
polynyas  and  then  drifted  with  little  change  in  thickness  toward  the  open 


ocean  (Pease,  1980; 

Hibler  and 

Ackley, 

1983)  . 

What  MIZ 

processes  are 

responsible  for  this 

des  true  tion 

and  how 

is  it 

coupled  to 

the  ocean  and 

a  tmosphere  ? 

Andreas  et  al.  (1984)  have  made  a  data  study  of  the  modification  of  the 
atmospheric  boundary  layer  in  the  MIZ.  During  an  episode  of  on-ice  winds  in 
the  Weddell  Sea,  they  obtained  wind  and  temperature  soundings  at  the  ice  edge 
and  four  additional  temperature  soundings  along  a  150  km,  along-wind  track 
into  the  ice  cover.  In  their  analysis,  they  have  suggested  that  the  destruc¬ 
tion  of  sea  ice  in  the  MIZ  is  episodic.  In  particular,  they  have  suggested 
that  during  episodes  of  strong  on-ice  winds,  there  is  a  large  air- to-surface 
heat  flux,  perhaps  due  to  secondary  flows,  which  could  result  in  rapid  ice 
abla  tion. 

A  two-dimensional  model  of  the  atmospheric  boundary  layer  based  on  the 
Herman  and  Goody  (1976)  arctic  stratus  cloud  model,  together  with  simple  sea 
ice  and  oceanic  mixed  layer  thermodynamics,  is  used  to  make  a  case  study  of 
the  Andreas  et  al.  (1984)  data.  We  use  their  ice  edge  wind  and  temperature 
soundings  as  a  prescribed  inflow  and  calculate  the  modification  of  the  atmo¬ 
spheric  boundary  layer  as  it  flows  on-ice.  A  model  simulation  which  incor¬ 


porates  the  exchange  coefficients  appropriate  for  a  rough  broken  ice  cover 
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(Andreas  et  al.,  1984)  is  used  to  test  several  results  and  hypotheses  of  the 
Andreas  et  al.  (1984)  analysis.  This  model  simulation  is  also  compared  to  a 
simulation  which  uses  prescribed,  open  ocean  exchange  coefficients  in  order  to 
obtain  an  estimate  of  the  effect  of  incorporating  realistic  surface  roughness 
on  calculations  of  the  surface  wind  stress  and  heat  flux.  These  simulations 
are  referred  to  as  the  "rough"  and  "smooth"  simulations,  respectively. 

Background 

The  Andreas  et  al.  (1984)  data  study  has  several  specific  results  and 
hypotheses  that  can  be  readily  tested  with  our  model.  The  principal  results 
and  hypotheses  include: 

a)  observed  soundings  and  height-fetch  cross-sections  of  potential 
temperature  that  can  be  compared  with  model  results; 

b)  there  is  a  large  air- to-surface  heat  flux  over  the  rough,  MIZ  ice 
cover,  perhaps  due  to  secondary  flow,  which  results  in  large 
rates  of  ice  ablation; 

c)  the  rise  with  fetch  of  the  isentropes  above  the  inversion  base  is 
the  result  of  adiabatic  lifting  due  to  upward  motion  forced  by 
low-level  convergence. 

One  goal  of  the  Andreas  et  al.  (1984)  study  was  to  measure  the  drag  coef¬ 
ficient  over  the  MIZ  ice.  Overland  (1985)  has  recently  reviewed  a  number  of 
field  measurements  of  the  drag  coefficient  over  sea  ice.  In  his  review,  he 
emphasizes  the  dependence  of  the  drag  coefficient  on  form  drag,  the  upwind 
surface  roughness,  the  height  of  the  atmospheric  inversion  and  the  atmospheric 
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stability.  The  relationship  between  the  drag  coefficients  and  the  exchange 
coefficients  for  heat  and  moisture  has  been  discussed  by  Walter  et  al.  (1984). 

An  important  aspect  of  this  study  is  the  modification  of  the  air  flow 
over  a  change  in  surface  roughness.  Much  of  the  earlier  field  work  on  this 
topic  is  limited  to  very  shallow  layers  and  short  fetches.  These  studies 
include  Bradley  (1968);  the  Riso  studies  (Petersen  et  al.,  1980  and  references 
cited  therein),  and  Ogawa  and  Ohara  (1985).  The  Andreas  et  al.  (1984)  study 
is  the  first  MIZ-scale  study  of  on-ice  air  flow.  Off-ice  flow  has  been 
examined  by  Overland  et  al.  (1983)  and  Reynolds  (1984),  while  the  analogous 
passage  of  cold  air  over  the  Great  Lakes  has  been  considered  by  Lenschow 
(1973)  and  Stage  and  Businger  (1981).  Other  relevant  studies  are  Rao  et  al. 
(1974)  and  Taylor  (1971),  who  have  done  some  modeling,  and  Hogstrom  and 
Smedraan-Hogs trom  (1985),  who  have  discussed  some  data  and  modeling  work. 
Additional  data  and  modeling  work  is  being  done  in  association  with  MIZEX,  the 
Marginal  Ice  Zone  Experiment  (see  Johannessen  et  al.,  1983). 

Among  published  models  of  the  atmospheric  boundary  layer  applicable  to 
the  MIZ,  the  one  most  relevant  to  our  study  is  the  one  used  at  the  Pacific 
Marine  Environmental  Laboratory  (Overland  et  al.,  1983;  Reynolds,  1984).  It 
is  a  one  layer  model  and  is  used  to  model  off-ice  air  flow.  The  one  layer 
approximation  is  useful  during  off-ice  situations  because  strong  surface 
heating  of  the  cold  air  over  the  relatively  warm  water  leads  to  convection  and 
a  uniform  sounding  of  equivalent  potential  temperature.  However,  the  situa¬ 
tion  observed  by  Andreas  et  al.  (1984)  is  stably  stratified.  Our  approach  to 
determining  the  turbulent  fluxes  and  stability  to  secondary  flow  is  to  make 
multi-layer  calculations  of  the  wind  shear  and  temperature  lapse  rate. 
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Recently,  Davidson  et  al.  (1984)  have  used  a  single  layer  model  of  the  marine 
atmospheric  boundary  layer  which  permits  a  constant  vetical  gradient  in  the 
mixed  layer  properties.  These  models  are  thus  distinctive  in  their  approach 
or  applicability. 

Model  Description 

A  multi-level,  linear,  Boussinesq  model  of  the  atmospheric  boundary  layer 
which  closely  follows  the  model  of  Herman  and  Goody  (1976)  is  used.  The 
principal  variables  are  0£,  the  equivalent  potential  temperature;  r,  the  total 
water  mixing  ratio  (the  sum  of  the  water  vapor  and  liquid  water  mixing 
ratios);  and  the  wind  velocity  components  of  u,  in  the  on-ice,  x  direction;  v 
in  the  cross-ice,  y  direction;  and  w,  in  the  upward  direction.  This 
atmosphere  is  described  by 
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Other  variables  include  K,  the  vertical  eddy  diffusivity;  the  radiative 

heating  term;  rL,  the  liquid  water  mixing  ratio;  and  f,  the  Coriolis 
parameter.  These  equations  have  been  linearized  about  a  base  state,  on-ice 
wind  UQ  of  13  ra  s~^. 

A  number  of  processes  are  included.  The  calculation  of  QRAD  and  K  and 
the  specification  of  Wp,  the  drop  fall  speed,  and  other  drop  variables  follow 
Herman  and  Goody  (1976).  Using  a  single,  constant  drop  size  in  the  radiation 
scheme  is  a  potentially  serious  limitation  (Tsay,  Jayaweera  and  Stamnes, 
1983),  but  we  believe  the  model  is  still  useful  in  identifying  the  basic 
physics  of  atmospheric  boundary  layer  modification  in  the  MIZ.  In  these 
calculations,  we  assume  we  are  in  the  northern  hemisphere  and  the  date  is 
April  21.  Pressure  p  is  calculated  by  the  Boussinesq  approximation 
(Schatzmann  and  Policastro,  1984). 

The  boundary  conditions  are  chosen  so  that  we  can  make  a  case  study  of 
the  Andreas  et  al.  (1984)  data.  Figure  1  is  a  schematic  diagram  of  the 
situation.  The  wind  and  potential  temperature  profiles  of  the  inflow  are 
prescribed  using  several  linear  segments  to  approximate  the  ice  edge  data  of 
Andreas  et  al.  (1984).  They  did  not  have  any  useful  data  on  the  moisture 
profile.  We  have,  therefore,  examined  the  Arctic  stratus  cloud  data  of  Tsay 
and  Jayaweera  (1984)  and  have  rather  arbitrarily  prescribed  a  "reasonable" 
inflow  moisture  profile.  Modeled  and  observed  inflow  profiles  are  compared  in 
Fig.  2.  Note  the  wind  shear  minimum  between  400-500  m.  We  will  consider  it 
below.  The  atmosphere  above  2  km  is  prescribed  from  the  60°N  July  U.S. 
Standard  Atmosphere  (Dubin  et  al. ,  1966).  We  have  chosen  to  prescribe  the 

turbulent  fluxes  at  the  top  of  the  model  domain  (2  km)  rather  than  prescribe 


6 


0g  and  r  or  use  a  zero  flux  condition  because  these  fluxes  are  more  consistent 
with  the  observed  wind  shear  and  vertical  gradient  of  potential  temperature 
near  2  km.  These  fluxes  are  -0.018  deg.  m  s  and  -6.8  x  10  m  s  , 
respectively.  A  positive  flux  is  upward.  The  initial  ice  concentrations  are 
based  on  the  observations  of  Andreas  et  al.  (1984). 

Bulk  aerodynamic  formulas  are  used  to  calculate  the  surface  turbulent 
fluxes  of  sensible  heat,  moisture  and  momentum.  The  open  ocean,  25  m  drag 
coefficient  is  set  to  equal  to  1.07  x  10  .  It  is  obtained  by  assuming 
neutral  stability,  a  constant  momentum  flux  between  10  ra  and  25  ra  and  a  10  m 
drag  coefficient  of  1.25  x  10-3  (from  Fig.  17.2  of  Charnock  (1981)  with  a  10  m 
wind  of  7.5  m  s-^).  In  the  "smooth"  simulation,  we  use  the  open  ocean  drag 
coefficient  at  all  grid  points.  In  the  "rough"  simulation,  we  accept  the 
Andreas  et  al.  (1984)  estimate  of  how  the  drag  coefficient  varied  with  fetch. 
Using  their  Eqn.  4,  we  set  the  ratio  of  the  drag  coefficient  over  ice  of 
concentration  C,  CD(C),  to  the  open  ocean  value,  C^0  ,  to  be 


CD(C) 


1+2.4  tanh  (2.5  C)  . 


The  applicability  of  their  Eqn.  14  is  discussed  below.  The  exchange  coef¬ 
ficient  for  heat,  C^,  was  not  estimated  by  Andreas  et  al.  (1984).  In  a  study 
over  the  Bering  Sea,  Walter  et  al.  (1984)  measured  the  ratio  C^/ to  be  0.20 
-  0.28  over  rough  sea  ice  which  had  a  drag  coefficient  of  =  3.0  ±  0.6  x 
10"3.  They  noted  that  their  measurement  of  C^/C^  was  consistent  with  earlier 
results  that  is  1.1  over  a  smooth  surface  but  decreases  rapidly  as  the 
surface  roughness  increases  beyond  a  certain  critical  value.  The  peak  Andreas 
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et  al.  (1984)  estimate  of  CD  =  4.0  x  10~^  is  larger  than  the  CD  measured  by 
Walter  et  al.  (1984)  and  suggests  that  C^/Cp  ought  to  be  small  in  our  case 
study.  Hence,  we  set  C^/C^  =  1.1  in  the  "smooth"  simulation  and  in  the 
"rough"  simulation,  which  has  a  wide  range  of  ice  concentration,  C^/C^  =  1-C. 
The  latter  relationship  is  ad  hoc  but  incorporates  the  gross  dependence  of 
Ch/Cd  on  surface  roughness.  The  exchange  coefficient  for  moisture,  Cg  is  set 
equal  to  (Walter  et  al.,  1984). 

The  surface  incorporates  a  simple  formulation  of  sea  ice  thermodynamics 
(ablation)  based  on  the  0  -  layer  model  of  Semtner  (1976)  and  a  heat  budget 
calculation  of  lateral  ice  ablation.  For  simplicity,  the  sea  ice  has  no  snow 
cover.  There  is  only  one  mixed  layer  temperature  in  a  grid  length  and  it  can 
not  exceed  the  freezing  point  of  sea  water  if  sea  ice  exists.  A  50  m  deep, 
motionless  oceanic  mixed  layer  is  used;  salinity  is  not  considered.  Further 
development  of  the  mixed  layer  part  of  the  model  is  planned. 

It  is  a  two-dimensional  model  with  forty  50  m  deep  layers  in  the  vertical 
and  seventeen  10  km  long  intervals  in  the  horizontal,  perpendicular  to  the  ice 
edge  direction.  Radiative  heating  and  convective  adjustment  are  calculated 
once  every  12  raintues.  At  time  intervals  of  2  minutes,  a  smoother  with  s^  = 
0.5  and  S£  =  -0.5  is  applied  to  each  vertical  profile  of  Qg  (Haltiner,  1971, 
pp.  270-274).  The  use  of  a  smoother  is  necessary  in  order  to  remove  some 
ripple-like  features  from  the  simulations.  An  Euler  forward  time  step  is  used 
with  the  radiative.  Upwind  differences  are  used  to  calculate  the  lateral  and 
vertical  advection  terms.  Time-stepping  of  these  and  the  other  terms  in  Eqns. 
(1)  to  (4)  use  Method  A  of  Young  (1968)  with  a  time  step  of  40  s.  This  is 
also  the  time  increment  for  the  surface  heat  budget  and  ice  ablation  calcula- 
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tions.  Although  these  time  steps  were  satisfactory  for  a  number  of  runs,  the 
rough  simulation  required  a  time  step  of  13.3  s  with  Method  A. 

Results 

The  vertical  soundings  and  fetch-height  sections  of  potential  temperature 
of  the  rough  simulation,  Figs.  3  and  4,  respectively,  generally  agree  well 
with  the  observations  of  Andreas  et  al.  (1984),  their  Figs.  2  and  6, 
respectively.  We  find  the  base  of  the  inversion  rises  from  a  height  of  525  m 
at  the  ice  edge  to  about  875-900  m  at  a  fetch  of  150  km.  It  also  crosses 
isentropes  from  its  initial  potential  temperature  of  4°C  to  a  final  value  of 
about  6.5-7.0°C.  Andreas  et  al.  (1984)  observed  the  outflowing  inversion  base 
to  be  at  a  height  of  1050  m  and  a  potential  temperature  of  6.4°C.  The  top  of 
the  inversion  is  smoothed  out  in  the  model  simulations.  The  modeled  contours 
of  potential  temperature  within  and  above  the  inversion  rise  with  fetch 
although  the  modeled  cooling  is  smaller  than  observed.  Below  a  height  of 
500  ra,  there  is  a  small  warming  with  increasing  fetch  in  the  model  atmosphere 
but  a  cooling  of  perhaps  1C  in  the  observations.  The  gross  features  of  the 
modeled  cloud,  shown  by  the  stippling  in  Fig.  4,  agree  with  the  observations 
of  Tsay  and  Jayaweera  (1984).  These  features  include  a  cloud  top  at  the  base 
of  the  inversion  and  at  longer  fetches  the  occurrence  of  the  largest  liquid 
water  mixing  ratios  near  the  top  of  the  cloud. 

Although  this  model  is  not  capable  of  simulating  secondary  flows,  we  hope 
to  be  able  to  say  whether  they  are  likely  to  exist.  Our  approach  is  to 
associate  the  model  wind  and  temperature  profiles  with  the  mean  flow  on  to 
which  the  horizontal  vortex  rolls  of  the  secondary  flow  are  added  and  test 
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whether  this  mean  flow  is  stable  to  secondary  circulations.  There  are  two 
major  energy  sources  for  rolls  --  convective  and  dynamic  instabilities  (Brown, 
1980).  Theoretically,  these  sources  can  act  independently.  However,  Muller 
et  al.  (1985)  have  noted  that  the  rolls  are  generally  observed  in  situations 
with  unstable  stratification  and  these  two  sources  usually  can  not  be  distin¬ 
guished.  They  have  suggested  that  inflection  point  instability  acts  as  a 
trigger  for  roll  development  which  is  enhanced  by  buoyancy  if  there  is 
unstable  stratification. 

If  the  atmosphere  were  dry,  it  would  be  stably  stratified  since  ©  *  0. 

3 

The  subscript  z  indicates  yy  .  The  moisture  profiles  must  be  considered 
before  we  can  say  whether  the  soundings  of  equivalent  potential  temperature 
are  stable.  We  can  say,  however,  that  if  the  500-m  thick  layer  between  the 
heights  of  375  m  and  875  m  were  to  have  (©E)z  =  0,  then  the  water  vapor  mixing 
ratio,  ry,  would  have  to  decrease  by  approximately  1.3  g/kg  over  that  depth. 
In  the  model,  however,  ry  becomes  increasingly  uniform  below  the  inversion 
base  as  the  fetch  increases.  Since  the  air  is  warmer  than  the  surface  in  the 
model  and  they  are  approximately  the  same  in  the  observations,  we  expect 
little  buoyancy  to  be  generated  at  the  surface.  We  conclude  that  convection 
is  unlikely. 

Following  Brown  (1980),  we  test  for  dynamic  instabilities.  The  Reynolds 
number  Re  is  on  the  order  of  600  or  larger,  which  is  in  the  inflection  point 
instability  regime  (LeMone,  1973).  We  consider  a  longitudinal  roll  axis  at  an 
angle  e  to  the  geostrophic  wind  and  the  wind  profile  in  a  plane  prependicular 


to  the  roll.  If  the  local  Richardson  number 
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T  +  g/c 

Rer  ©  (-?  2  P) 

T  z 


at  the  inflection  point  is  greater  than  0.25,  then  the  inflection  point 
instabilities  are  suppressed.  In  the  above,  g  is  the  acceleration  of  gravity, 
and  T  is  a  mean  temperature.  Applying  this  test  to  the  rough  simulation 
gives  the  results  shown  in  Table  1.  The  local  Richardson  number,  Ri^,  is  much 
greater  than  0.25  at  all  angles  £  and  fetches.  The  flow  is,  however,  less 
stable  with  increasing  fetch  because  of  a  steepening  of  the  wind  shear.  We 
conclude  that  secondary  flows,  are  unlikely  to  occur  during  the  period  of  the 
Andreas  et  al.  (1984)  observations. 

Above  the  inversion  base,  the  rise  of  the  isentropes  with  fetch  is  the 
result  of  vertical  advection.  The  dominant  terras  in  the  heat  budget  at  these 
heights  are  the  cooling  due  to  vertical  advection  and  the  warming  due  to  lat¬ 
eral  advection  (see  Fig.  5).  In  the  reference  frame  of  an  air  parcel  flowing 
downstream,  lateral  advection  is  zero  and  vertical  advection  dominates.  The 
modeled  vertical  velocity  is  1  -  3  cm  s-^,  which  agrees  well  with  the  Andreas 
et  al.  (1984)  estimate  of  2  cm  s-^.  The  model  processes  forcing  the  upward 
motion  are  a  slowing  of  the  zonal  wind  because  of  a  Coriolis-forced  turning  to 
the  right,  which  has  a  maximum  at  a  height  of  about  0.5  km  and  at  fetches 
greater  than  80  km,  and  a  low-level  convergence  forced  by  surface  friction. 

The  heat  budget  for  the  rate  of  change  of  potential  temperature  of  the 
air  column  at  a  fetch  of  150  km  in  the  rough  simulation  is  shown  in  Fig.  5. 
The  corresponding  potential  temperature  sounding  is  shown  in  Fig.  3.  The 
turbulent  flux  is  downward  at  all  levels  but  a  minimum  wind  shear  at  800  m 


causes  a  flux  minimum  at  that  level  and  the  turbulent  warming  near  875  m. 
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Below  the  height  of  the  minimum  flux,  turbulence  acts  to  make  the  vertical 
gradient  of  potential  temperature  uniform.  This  region  is  the  planetary 
boundary  layer  and  it  is  capped  by  the  inversion.  Above  the  wind  shear 
minimum,  the  wind  fields  of  the  rough  and  smooth  simulation  are  the  same. 
Strong  longwave  cooling  is  associated  with  the  cloud  top  at  950  m.  Cooling 
due  to  adiabatic  lifting  is  also  strong  near  and  above  the  cloud  top. 

The  cloud  top/inversion  base  rises  to  a  greater  height  in  the  rough 
simulation  than  in  the  smooth  simulation.  The  effect  of  turbulence  generated 
by  surface  friction  extends  to  a  greater  height  over  the  rougher  surface. 
Hence,  the  height  of  the  wind  shear  minimum  and  the  depth  of  the  boundary 
layer  increase. 

The  effect  of  incorporating  a  realistically  rough  surface  into  calcu¬ 
lations  of  the  surface  wind  stress  can  be  estimated  by  comparing  the  rough  and 
smooth  simulations.  The  25  m  bulk  exchange  coefficient  for  heat  and  momentum 
of  the  two  simulations  are  compared  in  Fig.  6.  In  the  rough  simulation,  the 
decrease  in  the  25  m  wind  speed  across  the  model  domain  (Fig.  7)  agrees  well 
with  the  observed  decrease  in  21  m  wind  speed  from  10  ra  s-'*'  at  the  inflow  to 
8.0  m  s_I  at  the  outflow.  When  the  drag  coefficient  is  fetch  dependent,  the 
maximum  wind  stress  does  not  necessarily  occur  at  the  point  of  maximum  wind 
speed  or  maximum  drag  coefficient.  In  the  rough  simulation,  the  maximum  wind 
stress  occurs  upstream  of  the  maximum  drag  coefficient  because  of  the  cumu¬ 
lative  slowing  of  the  wind  across  the  domain.  The  cumulative  slowing  of  the 
wind  is  much  smaller  in  the  smooth  simulation.  The  domain-averaged  wind 

_  O  _  T 

stress  is  0.208  J  m  in  the  rough  simulation  and  0.129  J  m  in  the  smooth 


Simula  tion. 


The  rough/smooth  ratio  of  the  domain-averaged  wind  stresses 
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(1.61)  is  reduced  by  the  relatively  small  wind  stresses  at  fetches  greater 
than  110  km  despite  the  large  drag  coefficients.  These  results  are  relevant 
to  the  problem  of  estimating  MIZ-scale  ice  drift. 

There  are  similar  effects  on  the  surface  turbulent  heat  flux.  The  down¬ 
ward  sensible  heat  flux  in  the  lowest  several  hundred  meters  is  larger  in  the 
rough  simulation  because  of  the  larger  wind  shear.  At  long  fetches,  the  sur¬ 
face  heat  flux  decreases  because  of  the  smaller  bulk  exchange  coefficients  and 
slower  wind  speeds  of  the  rough  simulation.  As  a  result,  the  temperature  of 
the  lower  boundary  layer  and,  hence,  the  air-surface  temperature  difference 
increases.  This  warming  is  not  found  in  the  observations,  which  limits  the 
relevance  of  our  heat  flux  calculations  to  estimating  what  would  be  observed. 
However,  we  can  note  that  maximum  surface  heat  flux  occurs  near  the  middle  of 
the  domain.  Comparing  the  surface  wind  speed,  bulk  exchange  coefficient  and 
air-surface  temperature  difference,  we  find  the  surface  heat  flux  at  a  fetch 
of  130  km  is  smaller  in  the  rough  simulation  because  of  the  slower  wind  speed. 
These  results  indicate  that  quantitative  estimates  of  the  surface  turbulent 
flux  of  sensible  heat  and  momentum  must  incorporate  not  only  an  appropriately 
rough  exchange  coefficient  but  the  effect  of  the  surface  roughness  on  the  sur¬ 
face  wind  and  air-surface  temperature  difference  as  well. 

Discussion 

Our  model  estimates  of  the  surface  turbulent  flux  of  sensible  heat  are  an 
order  of  magnitude  less  than  the  estimates  made  by  Andreas  et  al.  (1984)  using 
the  integral  method.  The  dominant  term  in  the  cooling  of  the  model  boundary 
layer  is  the  longwave  cooling  at  the  cloud  top  and  we  believe  it  is  the  major 
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part  of  the  heat  loss  between  the  soundings  used  in  their  method.  Apparently, 
Andreas  et  al.  (1984)  did  not  consider  this  term  and  instead  attributed  its 
effect  to  a  surface  flux  of  sensible  heat.  Cloud  top  long  wave  cooling  has  an 
important  role  in  the  theory  and  modeling  of  cloud  topped  planetary  boundary 
layers.  It  has  been  discussed,  for  example,  by  Herman  and  Goody  (1976),  who 
used  a  model  upon  which  ours  is  based,  Lilly  (1968),  and  Nieuwstadt  and 
Businger  (1984). 

We  doubt  the  Andreas  et  al.  (1984)  equation  relating  ice  concentration 
and  drag  coefficient  is  generally  applicable.  As  we  noted  above,  the  drag 
coefficient  is  a  function  of  form  drag  (Arya,  1973;  1975),  the  upwind  surface 
roughness  (Macklin,  1983),  the  height  of  the  atmospheric  inversion  and  the 
atmospheric  stability  (Overland,  1985).  In  order  to  estimate  the  drag 
coefficient  in  MIZ  modeling  studies,  a  composite  table  of  drag  coefficients 
for  various  ice  regimes,  air  temperatures  and  other  variables,  such  as  that 
compiled  by  Overland  (1985,  Table  6),  would  perhaps  be  more  appropriate. 

How  reasonable  are  the  drag  coefficients  estimated  by  Andreas  et  al. 
(1984)?  We  suggest  that  the  drag  coefficients  are  actually  smaller  than  their 
es  tima tes . 

Although  a  number  of  field  measurements  have  shown  that  the  MIZ  can  be 
substantially  larger  than  the  open  ocean  CD,  the  Andreas  et  al.  (1984) 
estimate  of  =  4.0  x  10  over  ice  of  concentration  0.8  is  larger  than  any 
of  the  drag  coefficients  included  in  the  composite  table  compiled  by  Overland 
(1985,  Table  6).  The  ice  regime  of  small  rafted  ice  floes  with  a  concentra¬ 
tion  of  0.8  -  0.9,  characteristic  of  the  inner  MIZ,  has  among  the  highest 


given  in  his  table.  These  coefficients  are  repeated  in  our  Table  2.  The  ob- 
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served  meteorological  regime  is  T  ~  0°C  and  Z.  >  400  m,  where  T  is  the  air 
temperture  and  Zi  is  the  inversion  height.  The  larger  coefficients  in  Table  2 
are  for  Ta  <  -5  C.  In  this  situation,  the  open  leads  and  thin  ice  areas  are  a 
source  of  buoyancy,  which  increases  CD  (Overland,  1985).  Since  Ta  ~  0  in  the 
Andreas  et  al.  (1984)  data,  we  expect  little  or  no  buoyancy  to  be  generated  in 
these  areas.  In  addition,  stable  stratification,  which  is  characteristic  of 
these  data,  is  associated  with  smaller  drag  coefficients.  Inspecting  compo¬ 
site  table  of  Cp  found  in  Overland  (1985,  Table  6),  one  might  speculate  that 
CD  should  be  in  the  range  of  2.2  -  2.6  x  10-^. 

The  sensitivity  of  the  model  simulation  to  the  values  of  the  bulk  ex¬ 
change  coefficients  is  tested  by  making  an  "intermediate"  surface  simulation 
and  comparing  it  to  the  other  cases.  Since  we  speculate  that  the  Andreas  et 
al.  (1984)  estimates  of  are  too  high,  we  multiply  the  (25  m)  ratio 
Cp(C)/Cp°  by  a  factor  of  0.67.  In  order  to  increase  the  surface  heat  flux,  we 
set  =  C£  =  CD.  The  25  m  is  shown  in  Fig.  6.  The  surface  wind  speed  and 
wind  stress  are  similar  to  those  of  the  rough  simulation  (Fig.  7).  There  is 
again  good  agreement  between  the  observed  and  modeled  dependence  of  wind  speed 
on  fetch.  A  height-fetch  section  of  potential  temperature  (Fig.  8)  shows  that 
the  potential  temperature  of  the  lowest  500  m  at  long  fetches  agrees  better 
with  the  observations  than  the  rough  simulation  does  but  the  1°C  contour  still 
does  not  appear.  In  fact,  the  1°C  contour  does  not  appear  even  if  we  use  CD 
of  the  rough  simulation  and  set  CH  ”  CD*  As  expected,  the  adiabatic  lifting 
of  the  isentropes  above  the  inversion  and  the  rise  of  the  inversion  with  fetch 
are  smaller  in  the  intermediate  simulation  than  in  the  rough  simulation. 
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TABLE  2 


Over 


Drag  Coefficients  (10^  C^)  Referred  to 
Small,  Rafted  Floes  of  Concentrations  0.8 


10  m  Anemometer  Level 
-  0.9  (from  Overland,  1985). 


T  ~  0C 

T 

<  -5C 

T 

<  -5C 

a 

a 

a 

Zi 

<  300  m 

Zi 

>  400  m 

2.6 

3.0 

3.7 

Ta  =  air  temperature,  =  inversion  height 
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